in memoriam Claude Itzykson 



LPTENS-95/24 



Almost Flat Planar Diagrams 



Vladimir A. Kazakov 
Matthias Staudacher t 
and 

Thomas Wynter ' 



Laboratoire de Physique Theorique de 
l'Ecole Normale Superieure * 

We continue our study of matrix models of dually weighted graphs. Among the attrac- 
tive features of these models is the possibility to interpolate between ensembles of regular 
and random two-dimensional lattices, relevant for the study of the crossover from two- 
dimensional flat space to two-dimensional quantum gravity. We further develop the for- 
malism of large N character expansions. In particular, a general method for determining 
the large N limit of a character is derived. This method, aside from being potentially 
useful for a far greater class of problems, allows us to exactly solve the matrix models of 
dually weighted graphs, reducing them to a well-posed Cauchy-Riemann problem. The 
power of the method is illustrated by explicitly solving a new model in which only positive 
curvature defects are permitted on the surface, an arbitrary amount of negative curvature 
being introduced at a single insertion. 



June 1995 

t This work is supported by funds provided by the European Community, Human Capital and 
Mobility Programme. 

* Unite Propre du Centre National de la Recherche Scientifique, associee a l'Ecole Normale 
Superieure et a PUniversite de Paris-Sud. 



1. Introduction 

Hermitian one matrix models were introduced and for the first time solved in the 
large iV limit in the seminal paper by Brezin, Itzykson, Parisi and Zuber [Q. These mod- 
els generate ensembles of planar, random graphs whose vertex coordination numbers are 
controlled by the matrix potential. By varying the potential, different classes of diagrams 
may be obtained, e.g. random square or random triangular lattices. However, despite this 
freedom, there is a class of physically important lattices that can not be generated by 
simply tuning the potential: regular, flat lattices with fixed coordination numbers of both 
vertices and faces. To attain them it is necessary to study planar graphs having coordi- 
nation number dependent weights for both the vertices and faces. It is straightforward to 
define modified hermitian matrix models producing such graphs, but they can no longer 
be treated with the methods of UJ. In fact, until very recently this class of models of 
dually weighted graphs seemed intractable. However, an important but little noticed ob- 
servation due to Itzykson and Di Francesco @ has made possible the explicit treatment 
of dually weighted graphs. The number of degrees of freedom of these models is crucially 
reduced by rewriting the model in the language of group theory. It should also be noted 
that this method, based on expanding the matrix model potential in Weyl characters, was 
already used presciently in a special case in another early paper by Itzykson and Zuber U . 
In a recent work Q we demonstrated that this new approach leads indeed to a problem 
amenable to mathematical analysis once the large N limit is taken. 

The physical importance of matrix models has been elucidated through a large body 
of work over the last ten years. In || || matrix models were first introduced to furnish 
a description of two-dimensional quantum gravity and non-critical bosonic strings and 
successfully used to calculate the critical properties of these theories. This approach is 
based on the representation of the sum over world-sheet metrics as a sum over dynamical 
triangulations as originally proposed in || || J7|. Studying the crossover from random, 
dynamical graphs to regular, static graphs, then, will correspond to suppressing the curva- 
ture fluctuations of the world-sheet metric and result in a flat two-dimensional metric. Our 
work, in conjunction with 0, should thus be seen as representing a first attempt towards 
establishing a connection between integrable two-dimensional models both coupled to and 
decoupled from quantum gravity. 

To be precise, let us consider general planar graphs and introduce a set of couplings 
ti, <2, --t*, ... , namely the weights of vertices with 1, 2, q, ... neighbours, and a dual set 
t±, tii ■■■tq, ■■■ , the weights of the dual vertices (or faces) with appropriate coordination 
numbers. The partition function of closed planar graphs G is defined to be 

z(*v)=e n t*/ v;t * vq 

G v*,v q £G 
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where v* , v q are the vertices with q neighbours on the original and dual graph, respectively, 
and #v* , j^v q are the numbers of such vertices in the given graph G. Choosing t* = t q = 
5 q ^4 the only allowed graphs are regular square lattices (see Fig. l.a). 








(b) positive curvature ( c ) nega tive curvature 



Fig. 1: Flat space and curvature defects 



However, it is easy to see that a regular square lattice cannot be of spherical (i.e. pla- 
nar) topology. Positive curvature defects have to be added in order to be able to close the 
surface. Considering for the moment only even couplings, we must therefore "turn on" 
couplings t 2 or t%, or both (see Fig. l.b). Exactly four such defects are needed to close 
the square lattice. Adding more defects, then requires balancing the total curvature by 
also adding negative curvature defects. The simplest examples for such negative defects, 
corresponding to the couplings t§ and tg, are shown in Fig. I.e. Allowing for an arbitrary 
number of positive and negative curvature defects we expect to generate random graphs 
which, at critical values of the couplings, corresponding to very big graphs dominating the 
sum in the partition function ( |1 . 1|) , allow us to reach a continuum limit lying in the uni- 
versality class of pure two-dimensional quantum gravity || Q. On the other hand, having 
"tuned away" the negative curvature couplings t q , t* with q > 4, no such continuum limit 
is possible. Then, only a small, finite number of positive curvature defects are allowed; 
this brings us back to the phase of essentially flat surfaces. The main physical motivation 
for studying the models of dually weighted graphs, then, is to understand the transition 
between these two very distinct phases. 

In the present paper we continue to develop powerful techniques which permit us to 
address this physical problem. Furthermore, we will present the full and explicit solution of 
a non-trivial problem: the case of flat, planar graphs with an arbitrary number of positive 
curvature defects and a single negative curvature defect (see Fig. 2(a)) adjusted to balance 
the total curvature. We call the resulting lattice surfaces "almost flat planar diagrams". 
A typical surface of this type is shown in Fig. 2(b). 



2 



This model illustrates a non-trivial example that can be solved by the method pre- 
sented in this paper. This model cannot currently be solved by standard matrix model 
techniques 






(a) 



(b) 



Fig. 2 (a) Negative curvature defect of angle (2 — q)ir and (b) a typical surface. 

It should be stressed that the methods we develop here are general and could have 
applications going beyond the problem under investigation. Given that the model of dually 
weighted graphs seemed entirely inaccessible even a short while ago, we regard the present 
approach to be an important step in extending current large N techniques. 

We will quickly recall in the next section some of the results of our previous paper [Q 
and precisely define the class of models we are studying. Then, in section 3, we demonstrate 
how to derive the large iV limit of group theoretical characters. The model of almost flat 
planar diagrams will be solved and interpreted in section 4. The full model capturing the 
transition from flat to random graphs will be briefly discussed in section 5. We demonstrate 
how to reformulate it as a well-posed Cauchy-Riemann problem. We conclude in section 
6 and present an outlook on how our approach might be put to further use in the near 
future. Technical details and additional illustrations are included in two appendices. 

2. Review of the character expansion method for matrix models of dually 
weighted graphs 

The partition function (|1 . 1|) for dually weighted graphs can be formulated as the 
following matrix model: 




(2.1) 
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with 

Vb(MA) = V- TrB k (MA) k . (2.2) 

k=i 

The matrices A and B are fixed, external matrices encoding the coupling constants through 

= 5* and t g = lTr^. (2.3) 

The model generalizes, for A ^ 1, the standard one matrix model first solved by Brezin, 
Itzykson, Parisi and Zuber |]J. It can no longer be solved by changing to eigenvalue 
variables; a reduction to iV variables is nevertheless possible. An expansion of the potential 
into a sum over invariant group characters allows all integrations to be performed and ( |2 . 1| ) 
to be reformulated as a statistical mechanics model in "Young-tableau weight space" . This 
reformulation should be called, after its discoverers, the "Itzykson-Di Francesco formula" 
and reads 

Z(tX)=c £ n X{H}(B) (2.4) 

Here c is a constant that we can drop, the weights {h e } are a set of N/2 even, increasing, 
non-negative integers while the weights {h°} are N/2 odd, increasing, positive integers, and 
the sum is taken over all such sets. The characters can be defined through two equivalent 
formulae. The first is the Weyl formula: 

det. h n (a ! l l ) 

where the are the eigenvalues of the matrix A and A(o) is the Vandermonde determinant. 
The second definition makes use of Schur polynomials, P n (9), defined by 

e ^ lZ % = J2z n P n (e) with e % = - Tr[A*], (2.6) 

n=0 % 

in terms of which the character is 

X{h} (A) = det (kl) {P hk+1 - l (e)). (2.7) 

It was demonstrated in Q how to take the large N limit of this expansion. In this 
limit, the weights -hhi condense to give a smooth, stationary distribution dh p(h), where 
p(h) is a probability density normalized to one. For technical reasons we restrict our 
attention to models in which the matrices A and B are such that traces of all odd powers 
of A and B are zero. This means that the our random surfaces are made from vertices and 
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faces with even coordination numbers only. As was discussed in jl, this ensures that the 
support of the density p(h) lies entirely on the real axis, and thus simplifies the solution of 
the problemB. The matrix A will satisfy this condition if we introduce an ^ x y matrix 
\fa in terms of which A and the character X{h}(A) are given by 



A 



a 





and X{h}i A ) = X{t^}{a)x { h^i } {a) sgn[JJ(fcf - h°)], (2.1 



1,3 



We now focus our attention on three intimately related models which capture the 
transition from flat to random graphs. 

I. Va(MA) = V \ Tx[A 2k ] (MA) 2k . 

k=i 2k 

II. V Aa (MA) = - a (MAf. (2.9) 

oo 1 

III. V A (MA 4 ) = V — Tr[A 2k ] {MA A ) 2k 

Here A± is defined to satisfy Tr[(A4) fc ] = N5kA and A is as defined in (|2.8|). The first 
model is self-dual, i.e. vertices and faces having the same coordination number have the 
same weights. The second and third models are dual to each other (the lattice of one 
corresponds to the dual lattice of the other) and are in turn related to model I by a simple 
line map. That is, we place the diagonal of a square belonging to model III (or alternatively 
a four- vertex belonging to II) onto each propagator of model I. Thus the vertices and face 
centres of model I become the vertices of model III (or alternatively the faces of II). We 
illustrate this in Fig. 3 belowl. 



1 We do not want to suggest that models with odd coordination numbers cannot be treated 
with our methods. 

2 Note that this line-map is only valid on the sphere. The corrections of I and III will 
thus be different. A careful analysis shows that the spherical free energy of model I is precisely 
twice the free energy of models II and III (since there are two ways of choosing the diagonal of 
a square in III, or alternatively two ways of splitting a four- vertex of model II). Note also that 
this non-trivial correspondence is predicted from our formalism, since we indeed obtain the same 
N = oo equations in all three cases. 
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Fig. 3: Graphical relationship between models I, II and III 



From this line map one can see that the expectation values in models I and III are also 
the same. More specifically 

(±Tr[{MA) 3 % = (^Tr[(MA 4 f k ]} in (2.10) 

Notice, however, that they are not equivalent to (jjTr[(MA) k ])u in model II. 

We can now return to the discussion of the large N limit and write the saddlepoint 
equation for these three models. Looking for the stationary point in ( |2.4j) , one finds from 
0], in all three cases, the following equation, valid on an interval [6, a] with < b < 1 < a: 

2F(h) + -[ dti p-j- = -In h. (2.11) 
Jo h-h! 

The solution requires, evidently, the knowledge of the large N limit of the variation of the 
characters in eq. (|2.8|) : 

F(h k ) = 2-£- In { A 2 /; } ; . (2.12) 
v ; dh% A(h e ) v ; 

The determination of F(h) is the subject of the next section. Let us also recall here the 
definition of the resolvent H(h): 

H(h) = £dh> ML. (2.13) 

In 0] we demonstrated, via a simple functional inversion, how to relate the results of the 
weight formalism to the resolvent W(P) = (^Tr P ^ M ) of the matrix model ( |2.1j ). In 
the model investigated in this paper, however, it is more natural to study the correlators 
(jjTr((MA) 2q )) . The results of the following section will provide a simple way to calculate 
such moments. 
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3. Large N limit of the character 

In the saddle point equation ( |2.11|) we introduced the function F(h) defined in eq. 
( |2.12|) as the derivative of the logarithm of a character. This function F(h) depends upon 
the moments of the matrix A, i.e. it contains all the information on the weights that one 
assigns to the faces of our discrete surfaces. In order to proceed with the solution of the 
saddle point equation, one would like to take the large N limit of (|2.12|) and express F(h) 



in terms of H(h) (which specifies the Young tableau) and the set of moments t q of the 
matrix A (the weights assigned to the faces). 

In f4|] a contour integral formula relating H(h), F(h) and the set of moments t q was 
derived. We recall here a single essential step of the derivation, which we will make use of 
shortly. We observed that 



Tr[a q ] = 22 „ 2 — where K = K + 2 <lh,ki (3-1) 



and the matrix a is the ^ x y matrix introduced in (|2.8|). For notational simplicity we 



omit an index k on h. In the large N limit ( |3.1|) was then reduced to a simple contour 
integral 

h q = - / i^W)) where t 2q = ^-Tr[a% (3.2) 
q J 2m N 

Note that the definition of F(h) (|2.12|) differs from that in the derivation in || since we 



are now restricting our attention to the case where only the even moments of the matrix 
A are non-zerJi. 



As it stands, formula (|3.2|) is of little direct use. It can however be dramatically 



simplified as we sketch out below. We introduce a function G(h) defined as 

G{h) = e H ^ +F{ - h \ (3.3) 

in terms of which ( |3.2|) becomes 

t2 Q = - <f ^ G(h)«. (3.4) 
q J 2-ki 



3 Indeed, it might be asked why we do not directly use formula (3.5) derived in section 3. of 
[Q. There, the contour integration relation was derived for the general case where both even and 
odd moments are non-zero. However, in the special case where we then set all odd moments to 
zero, e F ^ contains a cut overlapping with the cut of e H ^ h \ In this case defining the contour 
encircling the cut of e H ^ is ambiguous. We have therefore rederived the result for the reduced 
case of only even non-zero moments. The same note of caution applies to formula (3.8) of H]. 
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Changing integration variables from h to G we arrive at 



where h(G) is the inverse of the equation for G(h) given in (|3.3j ), and the contour in the 
complex G plane encircles the origin. We now assume that there are only a finite number 
of non-zero couplings t q . We obtain immediately the solution: 

Q , 

q = l 

Here ip{G) is an as yet unknown function, analytic in the vicinity of the origin, with 



V'(O) = 0. It is trivial to see that this satisfies (3J3). Note that, strictly speaking, we can 
not solve equation (3J3) for q = since ( |3.4|) is not defined there. The 1 on the l.h.s. of 
) comes from the normalization of the density p{h) (See appendix A). 
The unknown function if)(G) is not fixed by ( |3.4| ) and depends on the specific model 



being studied. We now give a very simple physical interpretation to this function. Let us 
return to the Schur polynomial definition of the character ( |2.7| ) . Differentiating eq. ( |2.6| ) 
with respect to 6i we see that 

<} P n {6) = P n - q {6) with e q = ^-t 2q . (3.7) 



dd q " v ' « 2 q 

This implies immediately that 

il^ ln (^¥>(«))=g4^) ^ K = h t -2 qKt . (3.8) 
From (|2.1j), ( |2.4| ) and ( |2.8| ), we see that the left hand side of this equation is equivalent to 



differentiating the logarithm of the original matrix integral (|2.1| ) with respect to t2 q - In 
terms of the dual to this matrix integral (in which the weights t 2 q assigned to the faces are 
now the weights of the vertices) this is equivalent to differentiating the coupling constants 



of the dual potential. So, denoting the dual matrix by M, the left hand side of eq. (|3.8|) is 



equivalent to the expectation value (Tr(MB) 2q ). Now, comparing the right hand side of 
( p.8|) to equations (|3.1|) and ( |3.2|) , we see that we have the following relation in the large 
N limit, 

(iTr(MS) 2 ^) = j ^ h(G) G~\ (3.9) 



8 



G(h) being denned by (|3~3|). It is now simple to follow identical arguments to those used 
to simplify ( [3^4]) to ( |3.6|) to arrive at 



Q t 00 1 

h ~ 1 = E nk + Djv 1 * (^ s ) 2<? ) G<? - ( 3 - 10 ) 

g=l g =l 

Given G(h), we have, after a functional inversion, the correlators of the dual model. 

To find G(h) we have to connect eq. ( |3.10|) with the saddle point equation (|2.11|) . 
From (plaft we obtain 

fl-(fe) = In — h — + H(h) with H(h)= I dh'-^- : (3.11) 

where the first term on the right is the contribution from the flat part of the density, i.e. the 
empty part of the Young tableau. The integral from b to a is the contribution from the 
"excited" part of the density, i.e. the non-empty part of the Young tableau. Noting, from 
the definition of G(h) (|3.3|) , that InG(h) = H(h) + F(h), we replace the integral of ( |3.11| ) 
by the contour integral 

where the contour encircles the [6, a] part of the cut of H(h). The discontinuity across 
this cut is precisely ±z7rp(/i). Note also that F(h) has - at least for some range of the 
couplings - no cut on the interval [6, a] . If we now change the variables of integration from 
h to G, as previously, and shrink the contour in the complex G plane catching poles on 
the way (see appendix A) , we arrive at the following simple relationship between equation 
O and H(h): 

{-1\(Q-i)h Q 

e*to = U± -i[G q (h). (3.13) 

^ q=l 

Some words of explanation are in order to clarify the meaning of this equation. Inverting 
eq. ( |3.6| ) leads to a multi-sheeted function G(h). The general picture is illustrated in Fig. 
4. One of the sheets is the physical sheet and has two cuts, one corresponding to e H<yh \ 
the other to e F ^; we label this sheet G\(h). The sheets G2(h), ... , Gq(K) are all the 
sheets connected to Gi(h) by the cut of e F ( h '; there are exactly Q of these sheets, where 
Q is the maximum inverse power of G in ( |3.10| ). 
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Fig. 4: Analytic structure of G(h) 



In appendix B several examples are presented to illustrate explicitly this general analytic 
structure. 

Equation Q3.13|) together with (|3.10| ) contains sufficient information to find the loga- 
rithmic derivative of the character. These two equations represent a well-defined Cauchy- 
Riemann problem for F(h) which can be explicitly solved. We will present the solution 
elsewhere. 



4. Almost fiat planar diagrams 

We now have all the tools necessary to reduce our model of dually weighted graphs to 
a well defined Cauchy-Riemann problem. In this section we will analyse the case in which 
only positive curvature defects are allowed on the surface, arbitrary amounts of negative 
curvature being introduced at a single point. This is done by studying the particular case 
t q = £2^,2 + £4^,4, which generates the flat patches (see Fig. 1(a)) and the positive defects 
(see Fig. 1(b)). The correlators ( |2.10|) then correspond to the insertion of a single defect of 
curvature (2 — k) ir (see also Fig. 2). They will be extracted using eq.( [3.10|) , after explicit 
calculation of the function h(G). 

From the analysis of the large N limit of the character in the previous section, we 
know that the product in ( 3.13 ) contains only two sheets G q {h) (see Fig. 4). We label the 



physical sheet G\{h) and the sheet connected to it by the cut of e F ^ h \ G^CO- Taking the 
logarithm of equation ( |3.13| ), we summarize the information extracted from the large N 
limit of the character by 

F 1 (h)+F 2 (h)+H(h) = -ln(-A), (4.1) 

£4 
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where \nGi(h) — Fi(h) + H{h). The two sheets G\{h) and are glued together by 

the square root cut coming from F{h). The combination F±(h) + ^(/i), evaluated on the 
cut of F(h), is twice the constant part of F(h) on the cut (the discontinuous part of F(h) 
is of opposite sign on Fi(h) and F 2 (h) and is therefore canceled). We thus have the two 
equations 

U (4.2) 

2F{h) + ${h) =-lnh, 

the first coming from the large iV limit of the character ( |3.13| ) and the second being 
the saddlepoint equation ( |2.11| ). These two equations tell us about the behaviour of the 



function 2F(h) + H(h) on the cuts of F(h) and H(h) respectively. We have introduced 
the notation ]f(h) to denote the real part on the cut of F(h), and similarly for JjH(h). The 
principal part integral in ( |2.11| ) is thus denoted in ( f4.2|) by lfl(h). 

Our object now is to reconstruct the analytic function 2F(h)+H(h) from its behaviour 
on its cuts. To do this we have to understand the complete structure of cuts. First we 
notice from ( |3.6|) that G(h) is non zero everywhere in the complex h plane except at infinity. 
The combination F(h) +H(h) thus has no logarithmic cut point except for the one which 
starts from h = b. This corresponds to the end of the flat part of the density p(h). We 
introduce two functions F(h) and H(h) defined by 

F(h) = F(h) -Inh and H(h) = H(h) + In — (4.3) 



in terms of which ( [4.2|) becomes 

2F(h) + H(h) = ln( -Uh - b)) 

I „ (4.4) 

2F(h) + $(h) = ln(h-b). 

These two equations define the behaviour of 2F(h) + H(h) on all of its cuts. By standard 
methods we now generate the full analytic function 2F(h) + H(h). We introduce three cut 
points, a, b and c whose values are fixed by boundary conditions (the points a and b define 
the cut of H(h) and c defines the starting point of the cut of F(h) which goes from c to 
— oo) and generate the full analytic function by performing the contour integral 



2F(h) + H(h) = y/(h-c)(h-b)(h-a) 



ds ln(s — b) 



+ 



Uc H 2ttz (h-s)y/(s-c)(s-b)(s-a) 
ds ln(— t 4 (s — &)) 

c F 2ri(h-s)^Ts-c)(s- b)(s -^J 



(4.5) 
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The contours Ch and Cp are illustrated in Fig. 5(a). The slanted zigzag line corresponds 
to the cut of ln(/i — b). Expanding the contours, catching poles on the way and using the 
fact that logarithmic cuts have a discontinuity of ±m, we arrive at 



2F(h) + H{h) = In ^ + y/(h-c)(h-b)(h-a) 

1 1 



ds 1 



+ 



(h - s) ^/( s - c )(s-b)(s-a) 
a ds MnU 
b (h-s) ^/( s - c )( s -b)(s-a) 



(4-6) 

Fig. 5(b) clarifies the sign convention for a/ \h — c)(h — b)(h — a) on the real axis above 
and below the cuts. Note that, for the cuts of 1/ yj (h — c){h — b) (h — a) the signs on the 
cuts are inverted compared to Fig. 5(b), i.e. +i <-> —i. The integrals in (|4.6|) are defined 
to be along the upper side. 




+i 

www 



C 



_vp b —i a 

(wW\j 



+ve 



—i 



+i 



(&) Contours for 



(b) Sign convention for yjh — c)(h — b)(h — a) 



Fig. 5: Contours and sign conventions for (|4.5| ) and (^ 



To fix the constants a, b and c, we expand (^1]) for large h and compare the resulting 
power series expansion to that obtained from inverting (|3.6|) : 

2F( h)+ H( h) = ^ + -L^ + (^). (4.7) 

The terms of C(^t^) depend on the as yet unknown function ip(G). Expanding ( fl6| ) for 
large h and comparing to (|4.7| ) we find the two boundary conditions 

t 4 = q = e _7r "^ and — - %= = ^-y/a — c, (4.8) 

with K and K' complete elliptic integrals of the first kind, defined in terms of their re- 
spective moduli k and k' = \Jl — k 2 through 
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The first condition fixes k and hence the ratio of the distances separating the cut points, 
and the second condition fixes a — c , i.e. the scale. The condition needed to fix the 
position of the cut points along the real axis is provided by the condition that the density 
must be normalized to one. 

We now perform the integrals in Q4.6| ) and, after using the first boundary condition 
and an identity between elliptic functions^, we obtain 



2F(h) + H(h) = -lnh- ^sn" 1 (J SlJL fc) ( 4 .10) 

K V a — b 

where sn _1 (,2, k) is the inverse Jacobi elliptic function. Using the saddle point equation 
2F(h) + Jfl(h) = — Inh and the fact that the resolvent for the Young tableau can be written 
as H(h) = JjH(h) =F i^p{h), we can immediately write down the expression for the density 
of Young tableau boxes as 

«=k'(S'»' (4 - n) 

The Jacobi elliptic function sn(z, k) is a generalisation of sin(^) with quarter period 
K. In fact, in the limit k — ► 0, which corresponds to £4 — > 0, the expression for the density 
becomes precisely (2/tt) sin -1 (a/ (a — h)/(a — b)). 

Integrating p(h) from b to a and equating the answer to 1 — b to ensure that the 
density is normalized to 1 (the flat portion from to b gives a contribution 6), gives the 
final boundary condition 

t 2 

a = l + ^-(K 2 - EK), (4.12) 

where E is the complete elliptic integral of the second kind. 

From the expression for the density we now generate the full Young tableau resolvent, 
H(h), in the standard way and obtain following expression: 



h r a 

H{h) =- + / dti 

h-b J b 



h-h' 



. , /- 1 la-h , ( 9 4^< sn 1 
\nh-—sn~ 1 \ - + 21n 



(4.13) 



K V a-b \ q i/4[( a _ c ^ a _ b ^ 1/4 e 4 (0). 



Using the above expression for H(h), the eq. Q4.10|) for 2F(h) + H(h) and the quasi- 
periodicity of theta functions, we can write the expression for G(h), 



1 / 7T _i I a — h mK' 

Ha I 1 



4 For this and many other relations between Jacobi's elliptic functions and theta functions 
useful for performing the calculations of this section see e.g. ||],||. 
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and its inverse 

q 3 / 2 l0 1 [6~ 1 (-GD 
h = a --Q2 

where the constant D is given by 



04(0) 



(4.15) 



ivq 5 / 4 2q 5 / 4 K 1 

To simplify ( f4.15| ) we have used the definition of the Jacobi elliptic function in terms of 
theta functions. In view of eq. (|3.10|) , we see that we have now explicitly calculated the 
generating function for the correlators for models I and III. 

We will now expand eq.( [4.15|) and read off the correlators as the coefficients of the 
positive powers of G. Notice that ( [4. 15| ) is a multivalued function since the function 6±(z) 
is periodic as z is varied in the real direction and quasiperiodic in the imaginary direction. 
We must thus choose the correct zero of the #4 (z) function about which to expand. The 
physical sheet corresponds to expanding about the zero z = i2L ^~- Using the definition of 
the Jacobi elliptic function sn(w) in terms of theta functions and shifting the arguments of 
the theta functions using their quasi-periodicity, we rewrite ( f!.15| ) as the pair of equations 

h - a+ U e "mf with2thes olutionof f ^ e "°M- (4 - 17) 

Expanding this for small G we find that the first three terms give (as expected from Q3.10D ) 
h = ^2 + q +1 + 0(G). Expanding three orders further, permits us (using (|3.10|) ) to read 
off the first three moments of model III (which are also the moments of model I): 

< 1 Tr [(MA 4 ) 2 ] ) m =^(l + / 8 ) 

( 1 Tr [(MA4) 4 ] ) IH =T J-3 (-8(1 + / 3 ) + 3/| - Af 2 f, + h) (4.18) 

( 1 Tr [(MA 4 ) 6 ] ) IH =^^(81 + 90/s - 30/| + 40/ 2 / 3 + 10/f - 10/ 4 - / 5 ). 
where for convenience we have defined 
h JM = ^ - EK) 



4 (O) 7T 

h =J^ = ^ ((3 - 2k 2 )K 2 - 6EK + 3£ 2 ) 
fs =f| =±((2-^-3^) 

h =°j^ = ™? ( (6 _ 6 e + k^)K 2 - 10(2 - e)EK + 15£ 2 ) . 



(4.19) 
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and have then expressed these derivatives as combinations of the complete elliptic integrals 
K, E and their modulus k. 

We can now give a simple physical interpretation of these moments. The first two 
are directly related to the free energy J-{t2,t^). The latter is defined as the sum over all 
possible surfaces with the topology of a sphere that can be constructed out of flat space 
and positive curvature defects. It is impossible to put a flat surface onto the sphere, so 
positive curvature defects are needed to close the surface. Since the defects in this model 
have a deficit angle of 7r it takes precisely four of them to close the surface into a sphere. 
The surfaces are in the form of a cylinder with both ends flattened. The four t 2 defects sit 
at the corners. Below we illustrate the free energy for model III: 




(4.20) 



Note that the flattened ends can have an angle of twist between them. The four t 2 defects 
correspond to vertices %■ Tr [(MA4) 2 ], and all other vertices (with four legs) correspond 
to the vertex \ Tr [(Mi 4 ) 4 ]. We see that the first two moments can be written in terms 
of the free energy J-{t 2 -,t^ as 



< 1 Tr [(MA 4 ) 2 ] ) IU =2-^-F(t 2 ,U) 



1 

N 



( ^ Tr [(MA 4 ) 4 ] ) IH 



We thus read off the free energy 



Using (14T9D , the identity ^ 



T{t 2 i t 4 



2K 2 kk' 2 a 



t 4 

T 2 



192q- 



~ A k Ht2M) - 



(1+/3 



(4.21) 



(4.22) 



along with standard identies for differentiating 
complete elliptic integrals with respect to the modulus k, it is trivial to verify that the 



moment ( i Tr [(MA 4 ) 4 } ) in given in (|4.18| ) is indeed four times the derivative of the 
free energy with respect to £4 = q. 

Using the definition of ^3 in terms of derivatives of the first theta function, 6i(z), 
along with the standard definition of the theta function as an infinite product, allows us 
to write the free energy as 



-4_d_ 

8 dq< 



In 



00 

[n< 

n=l 



1-q 



2n> 



E 




(4.23) 
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In this form we recognize the argument of the logarithm to be the partition function for 
the torus. The derivative operator acts to mark a single point. We have thus found, as 
illustrated in equation ( }4.23| ), that the free energy can be written as the free energy for a 
marked torus. Below we illustrate the connection between a marked torus and the flattened 
cylinder diagrammatically. 




Fig. 6 Diagrammatic connection between marked torus and flattened cylinder. 



Starting from the mark on the twisted torus, flatten the torus across its width (this defines 
a point on the opposite side) , then cut along the flattening and open out the crimped torus 
into a cylinder with flattened ends. The two points at either end of the flattening on the 
torus become the four t 2 defects of the flattened cylinder. 

Higher order moments correspond to inserting a single negative curvature defect. The 
lowest order defect of this type is the insertion of negative curvature of deficit angle — 7r 
introduced by the vertex Tr[(MA 4 ) 6 ] (see Fig. 1(c) and eq. ( |4.18|) ). 

As a final check of our solution, we expand ( -k Tr [(Mi 4 ) 6 ] )/// in powers of q 
(this can be done directly from the expression for the moment in terms of theta functions 

( 1 Tr [(MA 4 ) 6 } ) in = 4 ( 9q 2 + 27q 4 + 81q 6 + . . .). (4.24) 

It is easy to verify that this correctly counts the number of diagrams. 

Further moments can be calculated by expanding ( |4.17| ) to higher order. They can 
always be written as sums of products of complete elliptic integrals. 

We now look for a continuum limit in which the size of graphs tends to infinity. One 
can see that the critical point, at which the size of the graphs diverges, is at q = 1. Note, 
however, that (since the critical q is 1), in stark contrast to two-dimensional quantum 
gravity || ||, the leading behaviour for the growth of diagrams is not exponential but 
merely power-like. To extract the power, notice that the product n^°=i(l ~~ Q 2n ) in ( |4.23| ) 
can be written in terms of the r\ function as q ,_1,/12 ?7(zr) where q = e~ 7TT . Making use of the 
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modular invariance of the rj(r) function under the the modular transformation r — > 1/r, 
we extract „ . 

x 7T 2 tj 1 

where we have defined a "continuum cosmological constant" \x through q = 1 — fj, i.e. 
\i = 7tt -f- C(r 2 ). We thus see, changing from fixed cosmological constant to fixed area 
by Laplace transform, that the number of graphs grows as a linear power of the area. 
Employing the conventions of quantum gravity |J 0, this would formally correspond to 
a "string susceptibility" 7^ = 4. 

We can easily understand this result by performing the calculation directly in the con- 
tinuum limit. We thus integrate over cylinders of all possible lengths t and circumferences 
s weighted with a factor s (corresponding to the modular twist between the two flattened 
ends) and a delta function for the area so as to count the number of surfaces of a given 
area A. We thus perform the following integral: 

r°° f°° dt 

T~\ dt ds s 6(ts - A) = A (4.26) 
Jo Jo t 

We see immediately a linear dependence on the area A, but also a divergently large con- 
tribution coming from small t. The most important contribution comes from the cylinders 
which are infinitely short and thus have the maximum amount of entropy coming from the 
modular twist. 

It is interesting now to investigate the behaviour of the correlators in the large area 
(q — > 1) limit. Quite generally, for matrix models the correlators correspond to surfaces 
with a boundary of length proportional to the power of the correlator, and one seeks a 
continuum scaling limit for very long boundaries. The correlators ( -k Tr [(MAi) 2n ] ) in 
the present model, however, introduce point-like negative curvature and we cannot look for 
a scaling limit involving long boundaries!. Nevertheless we can find the limiting behaviour 
of these negative curvature insertions in the limit of large area. 

Using the modular transformation r — > 1/r (with q = e _7rr ) for the formula ( f4.15|) , we 
can also extract a scaling limit for the generating function for the moments. Specifically, 
for the theta functions 0±(z) and 9^(z), we find that the dominant contributions are 



6Az) =-^e-^V 1 / 4 (coshy + 0(q' 2 )) 



1 



e x {z) =-^e-^V 1/4 (sinhy + 0(q' 2 )) 



(4.27) 



5 In principle, it is possible to study boundaries of arbitrary length by taking correlators of 
( jr Tr [(MA|) 2n ] ), which correspond to a boundary in the form of the end of a cylinder. 
Technically, however, we do not at present have the means to calculate such quantities. 
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where y = z/r and q' = e~ ^ . Holding y fixed as a parameter of order 1, we take the limit 
as r — > (corresponding to 5 — > 1) and work to the first two orders in r. Remembering 
that q = 1 — iTT + . . ., the constants D and a to the first two orders in r are given by 

B = t 2 -^' 1/4 (l + ^r) and a=1 + ti [_L_ ( i__^)I], (4.28) 



r 



We can now define a natural rescaled parameter x = -^p- and perform the inversion of the 
theta function to the first two orders in r to find the generating function for the correlators. 

The lowest order term gives the contribution in ( |3.10| ) and also a part that cancels with 

iz. 

a 



a — 1. The next order gives the contribution -A along with the generating function which 



we read off as: 



E< > KMA^] ) G - = I p4^Zf + ^ - - j with .= ^.(4.29) 

n=l 



-1 ^. ^. /'^-l^ ^.2 



Gt 2 



This has a simple square root singularity at the point x = 1. The series expansions for 
sin -1 x and (sin _1 (x)) 2 then give us the dominant contribution to the correlators in the 
large area limit: 

{7r 2n (n!) 2 
<^2n+l — 2 4 ™+ 4 ((n+l)!) 2 (2r l +3) ' 

where we have again introduced the parameter \i defined by q = 1 — \i i.e. (i = ttt + 0{t 2 ). 
The number of surfaces of fixed area A for a correlator ( Tr [(MA4) 2n ] ) is thus seen 
to be of the order of A n , with entropy coming from modular twists analogous to those for 
the free energy. The rather curious structure of the considered surfaces is thus evident - 
they consist of cylindrical "fingers" growing out from the negative curvature defect (see 
Fig. 2(b)). The square root singularity at x = 1 means that there is a tree-like growth of 
the number of ways to attach the fingers to their base at the negative defect. As in the 
case of the free energy, modular integrations cause a filamentary structure of very long 
cylinders to dominate in the large area limit. 



5. The onset of quantum gravity: Adding negative curvature defects 

The introduction of arbitrary numbers of negative defects, specifically alongside 
the components £4, t\ and positive curvature defects, £2, gives us a model in which we 
can tune away the curvature fluctuations of two-dimensional quantum gravity. The large 
N limit of the character in section 3 allows us to understand the analytic structure of 
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the solution and thus reduce the model to a well defined Cauchy-Riemann problem. The 
function G(h) now consists of two sheets below the physical sheet (see Fig. 4). An extra 
sheet which we label G^(h) is now attached to the sheet G^/i) of the previous section by 
a square root cut. We thus have the following two equations 

2F(h)+fl(h) =-\nh 

F 1 (h) + F 2 (h) + F 3 (h) + 2H(h) = - In A ' ( °' ' ' 

The first is the saddle point equation ( |2.11| ). The second comes from the logarithm of 
eq. (|3.13|) , where we define, as before, Fi(h) by lnGj(/i) = H(h) + F{(h). Along with the 
boundary conditions provided by the coefficients of the negative powers of G in ( |3.10|) , the 
system of equations ( |5.1| ) completely determines the solution to this problem. 



6. Conclusions and outlook 

In the present work we have demonstrated that our technique of character expansions 
for large iV matrix models may be successfully applied to the study of a novel, up to now 
inaccessible phase of almost regular planar diagrams. This required determining - quite 
generally - the large N limit of Weyl characters through the functional equation ( |3.13| ). 
Specializing to almost flat graphs, we have then found the exact generating function ( |3.10|) , 
( |4.15| ) of planar square lattices endowed with a single negative curvature insertion balanced 
by a number of positive defects. 

We feel that our observations could trigger the investigation of many new phenomena 
in two-dimensional physics and the combinatorial theory of planar graphs. However, most 
urgent is the understanding of the crossover phenomenon from the phase of almost flat two- 
dimensional space to the phase of two-dimensional quantum gravity. It requires the careful 
analysis of the well-posed Cauchy-Riemann problem of the last section. This investigation 
is pending. Aside its obvious mathematical interest, the solution of this problem could 
help to solve the hitherto inaccessible problem of R 2 quantum gravity in two dimensions. 

In addition to the even lattices considered in this paper, our methods allow the study 
of the "melting" of more general regular, or almost regular, lattices; e.g. triangular lattices. 

It is well known that there are many intriguing relations between integrable two- 
dimensional models on regular lattices and dynamical planar random lattices. It is tempt- 
ing to try to unify the two classes of models, a project one might term GUTi- Our work 
should be considered a first attempt into this direction, even though it must be noted that 
further methods will have to be developed in order to successfully treat matter coupled to 
dually weighted graphs. 
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Some of the results presented above could be interpreted as insights into the structure 
of the group SU (oo) (see section 3 on the large iV limit of Weyl characters.). Further 
insights into this direction might prove very useful for the treatment of higher dimensional 
matrix models, e.g. the principal chiral field, discrete string theories in physical dimensions 
and, one hopes, QCD. 
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Appendix A. Derivation of the inversion formula 

We start by proving that the constant coefficient in ( |3.6| ) is equal to 1 (the normal- 
ization of the density). To correctly normalize the density, p(h), we have to ensure that 



Using the fact that InG(h) = H(h) + F(h), we replace the integral by the contour integral 



with the contour Ch encircling the [6, a] part of the cut of H(h), as shown in Fig. 7. The 
zig-zag line corresponds to the logarithmic cut starting at h = b. Note that this is not a 
closed contour since at b there is a discontinuity across the cut of ±in. Evaluating G(h) 
around this contour we see that its argument goes from +in at h = b (below the cut) all 
the way around to — m at h = b (above the cut). We now change integration variables 
from h to G, with (in light of the comment above) the contour Cg in the complex G plane 
encircling the origin (see Fig. 7): 



1 - b = - L i ft,(G) inG(ft) = - L £ ^ ( " (G) in G) - • (A - 3) 



where h(G) is defined through ( |3.6|) . The contour starts and finishes on either side of the 
cut generated by lnG illustrated in Fig. 7 by a zig-zag line. The total derivative term picks 
up the discontinuity across the cut giving —b. The final term, which picks up the constant 
coefficient of h(G), is thus equal to 1. 




(A.l) 




(A.2) 
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Fig. 7: Definition of contours Ch in the complex h plane and 
Cq and Cq in the complex C plane. 



We now complete the derivation of ( |3.13| ). As discussed in section 3. we start by 
generating H(h) (related to the full resolvent by H(h) = H(h) +ln ) from the contour 
integral 

1 ' J Ch 2ni h-h x 1 ' 

Changing integration variables from h to C, as above, this can be written as 

1 dG , „ ft'(G) 



where h(G) is defined through (|3.6|) and /i'(G) is the derivative with respect to G. 

We now simplify this contour integral by evaluating it for large h. Knowing the 
solution in any neighbourhood of h means that, by analytic continuation, we know it 
everywhere. For large enough h, we see from ( |3.6| ) that the contour in ( |A.5| ) will encircle 
precisely Q zeros of h(G), the zeros corresponding to the inverse powers of G. If we shrink 
the contour in ( |A.5| ) so that the contour hugs either side of the cut (see Fig. 7, contour 



Co) we pick up these Q poles: 

The remaining contour integral is relatively easy to evaluate provided careful attention 
is paid to the contribution coming from encircling the origin. The net result is that the 
contour Co contributes ln((— l)Q~ 1 t2Q) f rom encircling the origin and ln(h — b) from the 
discontinuity across the end points of the contour. Putting these results together and 
making use of the relationship between H(h) and H(h) we arrive at (|3.13|) . 

Appendix B. Analytic structure of G(h) 

As is discussed in Appendix A, the sheets G q (h) in the product of ( |3.13| ) are the 
physical sheet and all the sheets attached to the physical sheet by the cut of e H ^ h \ To 
clarify this we provide some simple examples. 
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B.l. Example 1. V B (MA) = 



In this simplest case it is immediate from eq.( |3.1(J| ) (since B = and thus i/j(G) = 0) 

that 

q=l 

This is a polynomial equation of degree Q. G(h) will thus be a multivalued analytic function 
with Q sheets. The different sheets are connected by square root cuts, represented in Fig. 
8 below by the vertical walls. 




Fig. 8: Analytic structure of G(h) for V B (MA) = 



Note that cube roots and higher order roots are just special cases of the above structure. 
For example, a cube root in the diagram above is generated when the two square root cut 
points touch. 

The G q (h) that enter the product in fl3.13| ) are precisely all the solutions, i.e. all the 
sheets. It then follows that 

e H{h) = j^, (B.2) 

which corresponds to a completely flat density p(h) = 1 with support [0, 1]. 

It is seen from equation ( |B.lj ) that at h = 1, G(h) becomes infinite on one of its sheets, 
so there is a pole at h = 1 on what we call the physical sheet. For ip(G) non zero, the 
positive powers of G "soften" this pole and stretch it into a cut. The cut corresponds to 
exciting boxes in the Young tableau. The next example illustrates this. 
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B.2. Example 2. V B (MA) = MA 

Here B = A\. It follows from eg. ( |3. 10| ) and a simple diagrammatic inspection that 



h-\ = Y^ + G 
^ Gi 

q=l 



(B.3) 



This increases the degree of the polynomial by one from the previous example, introducing 
an extra sheet. The pole that was at h = 1 has now opened into a cut (the cut of e H ( h ^) 
connected to this extra sheet (see Fig. 9). 



G, 




sheet of he H< ~ h ^ 
cut of e H ^ 
physical sheet 



Fig. 9: Analytic structure of G(h) for V B (MA) = MA 

The G q (h) that go into the product of eq.( |3.13| ) are the physical sheet and all the 
sheets below. We thus obtain 

Mh) = _JL_ (B.4) 



G-(hy 

where G~ (h) is the topmost sheet. 

B.3. Example 3. V B (MA) = (MA) 2 

By inspecting the moments of the dual model, we obtain 

Q 1 

q=l 

The sheet structure is still polynomial, but now, due to the symmetry G — ► G~ x of equation 
( p.5| ), the top sheets are the mirror image inverses of the bottom sheets. 
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Fig. 10: Analytic structure of G(h) for V B (MA) = {MA) 2 



Again, what was a pole at h=l has opened into a cut connecting the physical sheet 
to the mirror image inverses of the bottom sheets. 

The above three examples clarify the meaning of equations ( ft.lUp and ( pj.!3| ). A simple 
functional inversion developed in 0] allows us to relate H(h) to the resolvent, (Tr[ p _ M ]), 
of the matrix model. To verify the methods of section 3, we have directly calculated the 
matrix resolvent of these models using loop equations and simple diagrammatic arguments. 



24 



References 



[1] E. Brezin, C. Itzykson, G. Parisi & J.-B. Zuber, Commun. Math. Phys. 59 (1978), 35. 
[2] P. Di Francesco & C. Itzykson, Ann. Inst. Henri. Poincare Vol. 59, no. 2 (1993) 117. 
[3] C. Itzykson & J.-B. Zuber, J. Math. Phys. 21(3) (1980) 411. 

[4] V.A. Kazakov, M. Staudacher & T. Wynter, Ecole Normale preprint LPTENS-95/9, 

accepted for publication in Commun. Math. Phys. 
[5] F. David, Nucl. Phys. B257 (1985) 45. 
[6] V.A. Kazakov, Phys. Lett. B150 (1985) 282. 

[7] J. Frohlich, in: Lecture Notes in Physics, Vol. 216, Springer, Berlin, 1985; 
J. Ambj0rn, B. Durhuus and J. Frohlich, Nucl. Phys. B257[FS14](1985) 433. 

[8] P.F. Byrd & M.D. Friedman, "Handbook of Elliptic Integrals for Engineers and Physi- 
cists", Springer, Berlin, 1954. 

[9] D.F. Lawden, "Elliptic Functions and Applications", Springer, New York, 1989. 



25 



